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ABSTRACT 


The effect of oscillating flow on the turbulent intensity 
distribution in the turbulent boundary layer on a flat plate was 
investigated experimentally. 

A wind tunnel incorporating a set of rotating shutter blades 
downstream of the test section was used to create oscillating flow 
with amplitude variations from 6 to 110 per cent of mean velocity 

- | 
and frequencies ranging from 5 to 80 ore per second at a local 
Reynolds number of 1.5 x 10°. A vertical traverse of the turbulent 
boundary layer was conducted with a hot wire anemometer, and the 
mean velocity profile and the turbulent intensity distribution 
were recorded. 
Results indicate that there is no radical effect on the turbulent 


intensity distribution in the boundary layer due to frequency or 


amplitude changes in the oscillating freestream flow. 
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Ll. Introduction 


Although the unsteady turbulent boundary layer is a common 
phenomenon, very little work has been done to better the under- 
standing of this complex problem. The steady turbulent boundary 
layer is a non-linear mechanism, and when this situation is further 
complicated by a time-dependent mean freestream velocity, a full 
analytical solution is a formidable task. 

Karlsson (ref. 1) investigated the nonsteady turbulent boundary 
layer by separating the velocity fluctuations into periodic and 
random components. By averaging the oscillations over a complete 
period and substituting the results into the equations of motion, 
he obtained equations similar to the steady turbulent boundary layer 
equations with an extra fluctuating term. He predicted and experi- 
mentally verified that the oscillations had little effect on the 
mean velocity profile in the turbulent boundary layer. 

Klebanoff (ref. 2) has recorded the distribution of turbulent 
intensity in the turbulent boundary layer on a flat plate immersed 
in a steady flow of very low freestream turbulence. A similar study 
has also been conducted by Reichardt (ref. 9). 

Morrissey (ref. 5) has investigated the effects of large amplitude 
flow oscillations on the heat transfer rates from a flat plate. His 
results indicate no measurable influence on local heat transfer rates 
from the superposition of large amplitude fluctuations on the turbulent 
boundary layer. These results would suggest that the eddy diffusivity 
and thus the turbulent intensity distribution should exhibit similar 


behavior. 


Ja 


The objective of the present work is to investigate experi- 
mentally the effect of an oscillating mean stream on the turbulent 
intensity distribution in the fully established turbulent boundary 


layer. 
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2. Experimental Equipment and Procedure 
2.1 The Wind Tunnel 
2.1.1 General Description 
The oscillating flow wind tunnel in the Aeronautics 

Department Laboratories of the Naval Postgraduate School was used for 
this investigation. This is an open circuit tunnel with a test section 
24 inches square and 223 inches long. The tunnel is of metal and Lucite 
construction. The inlet section shown in Figure 1 has a contraction 
ratio of 16:1, which, coupled with three high solidity screens, results 
in very low freestream turbulent intensities. Measured steady free- 
stream turbulence levels in the test section varied from 0.25 to 0.37 
per cent of the freestream velocity. 

The wind tunnel is driven by two 100 horsepower Joy Axivane fans 
with adjustable blades. Located in front of each fan is a set of vari- 
able inlet vanes which serve as controls for the velocity in the test 
section. To reduce velocity the vanes pre-swirl the air in the direction 
of fan rotation, reducing the fan work input. The pitch of the fan blades 
was set to produce velocities of 88 to 250 feet per second. The mean 
test section velocity was measured with a pitot-static tube:and an 
inclined draft gauge. 

2.1.2 The Rotating Shutter Valve 
A shutter valve was used to superimpose an oscillating 
Sinusoidal component onto the mean freestream velocity. It consists 
of four rotating blades driven by a variable speed motor. Figure 2 
shows the shutter valve with 5-inch blades installed. The motor speed 
and pulley ratio selection determine the frequency of oscillation, and 


the blade width (and to some extent the frequency) determine the 
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amplitude. For the present work a range of frequencies from 5 to 
80 cycles per second and amplitudes from 6 to 110 per cent of the 
freestream velocity was produced by employing three blade widths 
and two pulley ratios. 

2.1.3 The Test Section 

The test section is constructed of 2-inch thick 
aluminum plates on top and bottom and 2-inch thick stress-relieved 
Lucite on the sides. One back panel of Lucite has been replaced 
with wood for convenience in mounting the model and instrumentation. 
Figure 3 is a photograph of the test section with the side panels 
raised. 
2.2 The Model 

The flat plate model was constructed of a slab of phenolic 
material 24 inches wide, 55 inches long and 5/8 inch thick. It was 
provided with a rounded leading edge in order to avoid leading edge 
separation. A 5-inch tapered spoiler which could be adjusted from 
outside the tunnel was used to trim the pressure gradient. The model 
and spoiler are shown in Figure 4. The traversing mechanism was 
attached to the underside of the plate, and the probe was extended 
up through the plate and bent forward and down to the upper surface 
at a point 27 inches aft of the leading edge. The traversing mechanism 
was designed to be mounted on the model rather than on the tunnel test 
section in order to avoid vibration problems produced by the oscillating 
stream. The traverse mechanism, shown in Figure 5, provided a resolution 


of 0.001 inch in displacement, measured with a dial indicator. 
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2.3 Instrumentation 
A transistorized, ten channel, constant temperature hot 

wire anemometer was employed in this work (Security Associates Model 
1000). This equipment is shown in Figure 6. The hot wire consisted 
of a 0.00015 inch diameter and 0.080 inch long tungsten filament. 
The hot wire circuit contains a built-in analog computer which produces 
an output which is Linear in velocity. Since the circuitry is D.C. 
coupled, the total instantaneous velocity is available at the output. 
Calibration of the apparatus was performed similar to the method used 
by Murphy (ref. 6). The turbulent intensity was measured with a 
Ballentine True RMS Voltmeter. An oscilloscope trace of the velocity 
wave form was recorded photographically. 

The frequency of oscillation was determined by a magnetic pickup 
mounted outboard of the upper shutter valve blade shaft. The output 
pulse frequency was read on a Berkeley decade counter. The instru- 
mentation system is shown pictorially in Figure 7 and schematically 
in Figure 8. 

2.4 Procedure 

After the plate was installed in the test section and adjusted 

to zero incidence, a static pressure survey was made to ensure the 
existence of a zero pressure gradient. Results of the pressure survey 
are shown in Figure 9. Next, conditions necessary for a fully established 
turbulent boundary layer were determined. A hot wire probe with tips 
bent down close to the surface well into the boundary was used to traverse 
the plate in a longitudinal direction. This device, shown in Figure 10, 
provided an indication of the transition point and the turbulent intensity 


along the plate at a constant vertical displacement of about 0.020 inch. 
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A test section velocity of 110 feet per second was found to cause 
transition 1/ inches from the leading edge and lead to a fully 
established turbulent boundary layer at the test point, 27 inches 
from the leading edge (Ry = 1.5 & ioe). The results of this 
preliminary investigation are summarized in Table lI. 

For each experimental run the boundary layer thickness was 
determined, and a vertical traverse of the boundary layer was then 
conducted, during which the mean velocity and the turbulent intensity 


were recorded. 
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TABLE I 


TRANSITION DATA 








Ue, x (in) RY me ion’ Regime 
100 10 Cee Z? Laminar 
he 059 Laminar 
16 - 846 Laminar 
18 1970 Near Transition 
I P02 Transition 
105 ike 5S Laminar 
15 icon Laminar 
18 1.000 Transition 
20 Me ZAG, Turbulent 
110 LS 875 Laminar 
16 933 Laminar 
Jey 1.000 Transition 
18 1,030 Turbulent 
Z7 L.o76 Turbulent 


(Test Position) 


Photographs of the oscilloscope traces were taken at various positions 
during the vertical traverse of the boundary layer. Oscillation 
frequencies of 5, 10, 20, 40, and 80 cycles per second were investigated 
for each of three sets of blade widths (2-inch, 3-inch, and 5-inch) 
corresponding to fluctuation amplitudes of 6 to 110 per cent of the 

mean freestream velocity. A summary of the experimental operating 
conditions and normalized frequency and amplitude parameters is presented 


in Table II. 
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TABLE IT 


EXPERIMENTAL OPERATING CONDITIONS 


Mean Velocity: 110 ft/sec Local Reynolds Number: 1.57 x 10° 
Run Width (in cps if a Intensity (% U 
1 2 5 Gao 0.07 2.0 
2 2 10 13-0 08 2) 
3 a 20 260) 07 iA 
4 2 40 D2 a0) . 06 220 
5 Zs 80 104.0 06 Lgegs) 
6 3 80 104.0 ELO 4.5 
7 5 40 20 Aiki 6.7 
8 3 20 LO) 24 6.3 
9 3 10 13ei0 16 4.6 
10 3 5 6.5 78 Dis) 
LD 5 5 6 i. 10 Ze) 
12 5 10 Lame P95 2 10 
LS 5 20 AN s60 20.5 
14 5 40 Ze .60 Lone 
15 > 80 104.0 0.60 LES: 
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3. Experimental Results and Discussion 


The experimental results are presented in Table III (Appendix). 
Figures 11, 12, and 13 summarize the effect of frequency on the 
turbulent intensity distribution in the oscillating turbulent 
boundary layer for each of the three amplitude ranges investigated. 
Also included in these graphs are typical velocity profiles for 
each of the amplitude ratios. The results are presented as the 
ratio of local intensity to freestream intensity (N.) because the 
freestream magnitude is clearly a determining factor on the intensity 
magnitude in the boundary layer. This normalized turbulent intensity 
appears relatively independent of frequency for each amplitude range. 

Curves for the mean intensity distributions for each of the 
amplitudes are compared with the steady flow results of Klebanoff 
(ref. 2) and Reichardt (ref. 9) in Figure 14. A check of steady flow 
results was in close agreement with those of Klebanoff. Although 
there is a slight trend toward decreasing turbulent intensity with 
increasing oscillation amplitude, the general shape of the curve is 
clearly preserved. In each case the pronounced maximum is in close 
proximity to the Wall. Moreover, there are no radical differences 
between the distribution results measured in oscillating flows and 
those measured in steady flow. One would assume such differences 
small when comparing flows with freestream turbulence intensity levels 
ranging from 0.20 to 25.0 per cent. Moreover, the data indicate that 
relative intensity is essentially unaltered very close to the surface, 
which is the primary region of importance in determining heat transfer 


rates. With some degree of certainty it may be concluded that frequency 
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FIGURE II 
DISTRIBUTION OF TURBULENT INTENSITIES 
N, = .07 R, =1.5x10° 
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FIGURE 12 : 
DISTRIBUTION OF TURBULENT INTENSITIES. 
Ne = .18 Rx =1.5x10° 


SZ 





y/§ 





0.5 1.0 1.5 2.0 2.5 3.0 


FIGURE |I3 


DISTRIBUTION OF TURBULENT INTENSITIES. 
No®.77 R, 1.5 x10§ 


33 


REICHARDT (Ref. 9) 
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KLEBANOFF (Ref.1) 
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FIGURE 14 
COMPARISON OF TURBULENT INTENSITY DISTRIBUTIONS 
VARIOUS FLUCTUATION AMPLITUDES AND 
STEADY FLOW RESULTS 
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and amplitude produce vem), little alteration mm the Cueeulent 
intensity distribution in the turbulent boundary layer. 

Figures 15, 16, and 17 are the oscilloscope trace photographs 
taken at four locations (y/§ = 1.0, 0.5, 0.02, and 0.005) in the 
boundary layer for each of the vertical traverses. An interesting 
phenomenon was observed in the boundary layer at Low frequencies 
(for example, Figure 15 for f = 5.0 cps). The effect of these weak 
oscillations is to produce layers in the turbulent boundary layer in 
which turbulent bursts are seen. This is indicative of a late 
transition regime, that is, not a fully established turbulent boundary 
layer. This conclusion is born out by the measured velocity profile 
(Figure 11) as well as the data reported by Miller and Fejer (ref. 4). 

Moreover, the hot wire traces are similar to those reported in 
the late transition regime by Miller and Fejer. This conclusion would 
also tend to explain part of the differences in results in Figure 14. 
The increased intensity observed for N, = 0.07 might be attributed 
to a late transition regime rather than a fully turbulent regime. 

Further work in this area might consider the effects on turbulent 
intensity profiles of a new velocity parameter involving the friction 
velocity. This approach might lead to stronger conclusions concerning 
the effects of freestream oscillations on turbulent intensity distri- 
bution. 

The RMS value of velocity measured in this work was a superposition 
of residual turbulence on the lower frequency oscillations. A more 
accurate method of obtaining true turbulence might be pursued in a later 
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4 » Cone lniswons 


The results of this investigation lead to the conclusion that, 
for the ranges of frequency and amplitude studied, the superposition 
of a fluctuating velocity component onto the mean flow does not 
substantially alter the turbulent intensity distribution in the 
turbulent boundary layer on a flat plate at zero incidence. This 
is in good agreement with the heat transfer measurements of Morrissey 
(ref. 5), who found only a 3 to 5 per cent increase in local Nusselt 
numbers in such flows. 

One concludes that the relatively unchanged turbulent intensity 
profiles are indicative of relatively unchanged eddy diffusivity 
distributions, confirming the analysis of Karlsson (ref. 1), who found 
that the mean properties of the turbulent boundary layer were unaltered 


by the imposition of freestream oscillations. 
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OL5 ~645 . 140 2.56 
Ok .600 aie 2.00 
.009 .10 mile 28 
7005 -400 . 168 3.08 
.003 -290 - 180 3 a8 
0.002 0. 130 02130 2 e338 


om =o om . ~ Fi ~—— ee eo)  cepenes siiseae SS SO NS SS SE th ES A  E 


PO IT A A OR ES aia A A Oe ee 
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TABLE III (continued) 





Run Number 7 N, = 0.16 f = 40.0 cps Tes =] 7, 
y/§ U/U eg u N. 
1.0 G77 70 0.070 1.04 

oie 985 078 1.16 
667 5 ie, 088 ie32 
. 200 -940 . L04 1.70 
416 evs: Ay es 4 1.96 
poe 900 ~toz 2628 
(220 . 880 - 166 2.48 
. 167 840 AO 2.64 
083 825 . 186 210 
050 yee, - 200 3.00 
,033 2739 2220 3.36 
O17 690 22 3.76 
ONL) - 665 ,298 3.86 
.O010 620 268 4.00 
007 “a6, 274 4.04 
.005 7/5 222 3252 
.003 3299 134 2200 
0.001 0.205 0.095 1.41 


a ag ES A A A EC ty NN et ee tT Sh RT RS et 
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TABLE III (continued) 








Run Number 8 N, = 0.24 f = 20.0 cps Te. = 6.3% 
y/§ u/U .. i N, 
Lo 02230 OnD7 3 Laghk6 

.833 .980 O81 128 
667 . 960 096 L52 
. 500 950 20 Lye 0 
:333 s925 ~159 22 
2250 895 m7 4 2a 
SUG? .875 . 176 20 
083 i350 . 180 ae 
050 .780 eo 3.04 
10382 740 201 3.34 
O27 .700 216 3.44 
.020 .670 4226 3460 
O17 615 423'0 3.66 
083 F570, 240 37180 
RO 118. 490 ~254 by SV 
. 007 425 2 6 C10 
.005 Fe a5. = lets) 3700 
0.003 Ou O39 220 


Rm mre Ait eet RP A A om RO nm ett 9 ED Dee Me ARN SN GR ES EG CRATE DET 2 nat CA A Ee NR tA A RY A ey RN i pa ont 
ee aes a eS me wet ae ae Oy ee eee ee 
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TABLE III (continued) 








Run Number 9 NW = 0.16 = 10 Of cns I. = 4.6% 
y/§ u/U . 1 N, 
UE, OR 70 0.056 P22 

1033 so7> 060 1.30 
667 . 960 076 1.65 
500 ps) 10) .096 2.08 
9333 900 122 2.64 
7250 880 by 2.7/6 
alk7 1a50 wli51 2.84 
.083 STAG ai32 2.88 
050 wf25 ob34 Zoe 
2033 1075 nat 3.06 
nOzs . 630 4150 3.26 
.020 ~ 580 . 168 3.66 
OLs 540 . 188 4.08 
1003 490 . 200 4.36 
~O010 410 - 180 3.90 
.007 320 E32 2.86 
005 2230 eal SLO 2.40 
0.003 0.210 0.104 2.06 
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TABLE [III (continued) 


Run Number 10 N. = 0 
y/& u/U ¢. 
0 One 

1055 »9e> 
Bicyey/ Ree) 
200 ne 0 
338 3920 
250 .880 
67 845 
083 .760 
067 ae) 
050 8016/7 
. 085 .625 
HOZ/, 200 
020 . 2210 
Ob 470 
.007 .300 
0.003 Omgb7 


20 


a en er me re ee Re me ei EM ne eR rg tn TE ay aE aN A Ms a 
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“A nt 


064 


.088 


110 


128 


. 0 


~ 144 


fe) 


ole 


BL, 


. 200 


7 oes 


. 204 


0 








TABLE III (continued) 





Run Number 11 N, = 1.10 t= 5.0geps oe = 25,07 
Ves u/U ,. I N, 
Ta) 05990 270 L508 

. 833 250 soo Pe 
.667 07 0 Cael, A s. 
oo 945 0 Le20 
416 5 els: .340 37 
5935 890 410 or 
ENG, .780 425 E78 
. 167 5750 435 1.74 
33 675 2450 L728 
. LOO . 630 462 1.85 
ye 0 491 Ley, 
5033 485 50) 20 
ZOLy 405 wooU 2 ol 
SOS 0 470 ee ei 
POLO nee 438 Lake 
.007 0 2518 1.41 
0.003 Oe00 Cr0> beee 


rts A: EE RT OP Sart, Aran RI RE a me arte eC Re i pte Se ey ily A eT OE 
ee ee eam ami aE rm mm A 
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TABLE III (continued) 


Run Number 12 N, = 0,95 f = 10.0 cps I. = 24.0% 
y/§ u/U,. I N. 
1.0 0.990 OF270 Lai 

, O50 .990 soe ee ks} 
667 77a ~285 1.18 
.500 . 965 pole) L230 
-416 .950 . 340 ea 
RS i 905 SoD Nou 
200) : 870 -400 1.67 
key .830 UG Level 
Bosh: . 800 red |S) 
FOO 120 -430 1.80 
. 067 Bo 15 5 Sis, V0 
Oss 960 300 2H 
, ON 4/0 . 940 2222 
Ols Be Ba 23015 
.010 . 300 ov U LyoT 
- 007 21 soo 0 2.28 
0.003 0.200 0.420 ir 


Ee RL SR re rm ut SRL OTE Oe ne OE ar at se rc NE Pe SE ee Ree gy RC at ceptng AN ME EO Nay A RT 
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TABLE III (continued) : 














Run Number 13 N, = 0.60 f = 20.0 cps Is = 20.5% 
we u/U,. I N, 
1.0 0.990 0.214 1.04 

eOOO 2700 PO ai ise: 1.06 
. 667 970 S229 1.10 
.500 - 940 A455) 135 
416 .920 .280 1338 
Pas foe .865 290 1.42 
250 805 . 300 1.46 
. 167 27.50 2320 <5 
e332 0 waco 1.60 
. LOO - 680 345 1.68 
.068 620 Po 8 Looe 
-033 280 ~420 2.05 
OFF -420 460 2500 
.013 £530 -440 2A 
.010 250 ~410 2.00 
007 AALS, -400 lee 
0.003 Oli LG, 0.360 La 
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TABLE III (continued) 








Run Number 14 N, = 0.60 f£i= 40.0 eps = 18.87 
y/§ u/U ,. if Fi 
EEO 0.990 0.225 20 

O28 wo5 BAA, S21 
.667 . 980 $235 a5 
.500 .960 .260 6 
416 Be es 276 47 
£333 875 2305 .63 
5 28)(0) .810 -330 ./76 
. 167 = TAG -340 iSd 
es .760 5 ENS) poo 
. LOO 2700 2355 .88 
.067 -ooW 270 7) 
.033 .600 425 226 
Oy 240 O50 .90 
013 -480 5500 75 
010 -400 oO A) 
007 .320 ae Ge -40 
0.003 0.200 0.216 5 
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TABLE III (continued) 





Run Number 15 N. = 0.60 £ = 80 cps I — 19.82 
vi S u/U.. it N. 
1.0 07.990 0.230 ns 

. 833 985 2235 Las 
667 A920 245 L423 
200 - 9510 260 aad 
-416 Fe) 25) 212 eis: 
p65 15: . 880 295 1.48 
290 810 .o20 1.62 
Gy ./60 - 340 es Il 
LoS 745 345 1.74 
. LOO 650 290 Lia 
.067 .615 «362 1,83 
2083 -490 .395 1.99 
O17 -440 430 Zale 
01s 420 415 ZG 
.010 385 . 365 1.85 
BOOT 20 347 Laws 
0.003 Ovid 0.297 1.50 
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